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Abstract: Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) fil- 
tering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a 
comparative study of different approaches to reduce this loss of accuracy within the framework of 
the dynamic Germano SGS model. One of the possible approaches is to apply Harten’s subcell 
resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid 
points adjacent to the exact shock location. The other considered approach is local disabling of the 
SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence 
interaction problem for comparison of the considered modifications of the SGS filtering procedure. 
For the considered test case both approaches show a similar improvement in the accuracy near the 
shock. 

Keywords: High order low dissipative scheme, Germano model, LES filtering, DNS, Subcell 
resolution scheme, Shock/Turbulence interaction 


1 Introduction 

The presence of a shock wave in turbulent flows might pose a numerical accuracy problem in employing 
the SGS filtered equations across shocks, depending on the LES model, the grid size, as well as the shock 
strength. Since the majority of dynamic LES models involve filter operations, hereafter referred to as “LES 
filters” to distinguish them from standard “numerical filters”, when the LES filtered equations are applied 
through the shock, the Rankine-Hugoniot relations are affected by the filtering operation, since the filtered 
variables are not discontinuous. In the present study we consider LES with implicit filtering employing the 
dynamic Germano procedure [1] for calculating the model coefficients. The dynamic Germano procedure was 
developed for shock free turbulence. Sagaut and Germano [2] have noticed that the usual filtering procedures, 
based on a central spatial filter that provides information from both sides, when applied around the shock, 
produce parasitic contributions that affect the filtered quantities. They suggested using non-centered filters 
to avoid this nonphysical effect. In [3] shock-confining filters have been proposed instead. Another approach 
based on the deconvolution method is considered in [4j. 

Aside from the subgrid scale filtering procedure, the accuracy of LES with shocks depends heavily on the 
accuracy of the numerical scheme. In this study we consider a combination of the low dissipative high order 
nonlinear filter scheme of Yee and Sjogreen [5] to locate the shock accurately, and the subcell resolution 
method of Harten [6] to confine the shock location to be within a grid cell. Previous studies indicate that 
the combination of the nonlinear filter scheme with Harten’s subcell resolution method is able to accurately 
locate the shock within a grid cell. For the comparison we will also consider a modification of the SGS 
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filtering procedure including local one-sided filtering without subcell resolution. One more modification 
considered here is local disabling of subgrid-scale (SGS) dissipation, which has been employed in previous 
studies [7, 8]. A new method to handle the transition points (buffer zone) between the smeared shock and 
the one-sided SGS filter equations is under development. 

The outline of the paper is as follows. Section 2 provides a description of the LES model considered and 
the numerical scheme used to solve the governing equations with the details on the three modifications of the 
LES filtering procedure. Section 3 describes the shock-turbulence interaction problem setup and a comparison 
of the the current DNS computations with the digitized results from [9]. Then LES computations employing 
the proposed method are compared with the filtered DNS data, and with LES results using standard and 
modified LES filtering procedures. 


2 Mathematical Formulation and Numerical Methods 

2.1 Governing Equations and LES Model 

We consider the filtered compressible Navier-Stokes equations written in the conservative form 

dtp + dj(pUj) = 0 (1) 

dt(pui) + dj (pui.Uj + pSij - f ij + rfj) = 0 (2) 

d t {pE) + dj(pEuj +pUj - fijPi +qj + qf) = 0, (3) 

where p is the density, u* is the i th velocity component, p is the pressure, T is the temperature, E is the 

total energy, and t is the time. For a function /, the LES filtering operation is denoted as /: 

i r°° r°° , x 

= x / / G(z£,t-T)f(Z,T)dSdT, (4) 

^ J — oo J — oo ' ' 


where G(x,t ) is the filter kernel in physical space and A is the filter width. The Favre filtering operation is 
denoted as / = pf/p, and / stands for the Favre-filtered variables: 


Tij — 2p,(T)(Sij fiij&k U'k) i ^ij — (Ay ll t “h 0; V. j ) / 2 


(5) 


Qj = -A (T)djf. (6) 

Unlike the “bar” and “tilde”, the “breve” symbol does not denote a filtering operation, but indicates that the 
quantity is based on primitive filtered variables. In the equations (5) and (6) the dynamic viscosity is given 
by p(T) = pq(T /Tq ) 3 / 4 and thermal conductivity is expressed through a constant Prandtl number Pr and 
heat capacity at constant pressure c p : X (T) = c p p(T) / Pr. The equation of state is p = RpT, where R is 
the gas specific constant. The subgrid-scale (SGS) terms, SGS stress tensor rA , and SGS heat flux qA are 
modeled as follows: 


T ij - o T kJij 


^~kk 

= 


2/ '/ ( S'i j ,, ‘Xk k G y ) 

2CipA 2 \S\ 2 

dj f , 

Pr* 3 


(7) 

(8) 

(9) 


where p* = pC s A 2 \S\, and |5| = \j2SijSij. The Smagorinsky constant C s and the constant for the isotropic 
part of the SGS stress Cj are obtained through Germano-Lilly [10] procedure: 


( L ij M ij s ) H c = (. Lu) h 

1 (K 1 ) ’ 
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where 



— ? ^ij — ^ P'U'iP'Uj / P 

(11) 


= — 2pA 2 2 (X- - + 2A 2 (p|S|4) - * (pISISuSij) 

(12) 

K 1 

= 2^A 2 |I| 2 -2A 2 (p|5| 2 ) 

(13) 


and < f >h stands for averaging in homogeneous (periodic) directions. 

For the considered test case with low turbulent Mach number M t < 0.4 it is shown [11] that the isotropic 
part of the SGS stress can be neglected: Cj = 0. The Germano procedure requires an explicit test filtering 
operation, denoted here with the hat symbol. For the Germano procedure we use a 3D filtering operator 
based on a ID trapezoidal filter: 

fi = J fi-l + 2 fi + j/<H fi- (14) 

Note that according to [12] the width of the discrete filter can be estimated as 

(N-l)/2 

A = ft 12 3 2W h ( 15 ) 

\ j=-(JV-l)/2 

where ft is the grid spacing and Wj are the filter weights. Using formula (15) for filter (14) gives A = h\J 6. 

2.2 High-Order Filter Schemes 

In order to solve the system (1) - (3) introducing as little numerical dissipation as possible, we use the 
high-order nonlinear filter scheme of Yee et al. [5, 13, 14, 15] which consists of three steps. 

2.2.1 Preprocessing Step 

Before the application of a high-order non-dissipative spatial base scheme, a preprocessing step is employed to 
improve the stability. The inviscicl flux derivatives of the governing equations are split into the following three 
ways, depending on the flow types and the desire for rigorous mathematical analysis or physical argument. 

• Entropy splitting of [16] and [17, 18]. The resulting form is non-conservative and the derivation is 
based on entropy norm stability with boundary closure for the initial value boundary problem. 

• The system form of the [19] splitting. This is a conservative splitting and the derivation is based on 
physical arguments. 

• Tadmor entropy conservation formulation for systems [20]. The derivation is based on mathematical 
analysis. It is a generalization of Tadmor’s entropy formulation to systems and has not been fully 
tested on complex flows. 

For the current test case containing a shock wave in the flow field it is more appropriate to use a conservative 
splitting. The Ducros et al. splitting is employed for all the computations. 

2.2.2 Base Scheme Step 

A full time step is advanced using a high-order non-dissipative (or very low dissipation) spatially central 
scheme on the split form of the governing partial differential equations (PDEs). A summation-by-parts 
(SBP) boundary operator [21, 22] and matching order conservative high-order free stream metric evaluation 
for curvilinear grids [23] are used. High-order temporal discretization such as the third-order or fourth-order 
Runge-Kutta (RK3 or RK4) method is used. It is remarked that other temporal discretizations can be used 
for the base scheme step. 
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2.2.3 Post-Processing (Nonlinear Filter Step) 

To further improve nonlinear stability from the non-dissipative spatial base scheme, after the application 
of a non-dissipative high-order spatial base scheme on the split form of the governing equation(s), the 
post-processing step is used to nonlinearly filter the solution by a dissipative portion of a high-order shock- 
capturing scheme with a local flow sensor. The flow sensor provides locations and amounts of built-in 
shock-capturing dissipation that can be further reduced or eliminated. At each grid point a local flow sensor 
is employed to analyze the regularity of the computed flow data. Only the discontinuity locations would 
receive the full amount of shock-capturing dissipation. In smooth regions no shock-capturing dissipation 
would be added. In regions with strong turbulence, if needed, a small fraction of the shock-capturing 
dissipation would be added to improve stability. For a variety of local flow sensors with automatic selection 
of the proper parameter, depending on different flow types see [13]. For the problem considered in this work 
we use Ducros et al. flow sensor [24]: 


(V-u) 2 

(V • u) 2 + w 2 + e 


(16) 


Here u is velocity vetor, u> is vorticity magnitude and e = 10~ 6 is a small number to avoid division by zero. 
The nonlinear dissipative portion of a high-resolution shock-capturing scheme can be any shock-capturing 
scheme. For the problem considered in this study it is activated when the Ducros et al. sensor w > 0.6 for 
the case M = 1.5 and w > 0.3 for the case M = 3. 

Let U* be the solution after the completion of the base scheme step. The final update of the solution after 
the filter step is (with the numerical fluxes in the y- and z-directions suppressed as well as their corresponding 
y- and z-direction indices on the x inviscid flux suppressed) 


U 


n+1 

j,k,l 


= u* 


At \H 

J.k.l - Xx [H : 


3 + 1/2 H j- 1 / 2 ] 


Hj- (-1/2 — Rj+1/2H j+1/2, 


(17) 


where Rj + 1/2 is the matrix of right eigenvectors of the Jacobian of the inviscid flux vector in terms of Roe’s 
average states based on U* . H* +1 , 2 and H*_ 1 , 2 are “filter” numerical fluxes in terms of Roe’s average states 

based on U*. Denote the elements of the filter numerical flux vector iJ J+1 / 2 by lij +i / 2 - 1 = 1,2,..., 5. The 

nonlinear portion of the filter numerical flux li 1+1 j 2 has the form 

hj+ 1/2 = 2 u J+ 1 / 2 < ^i+ 1 / 2 ' ( 1 ®) 

Here tol +1 / 2 is a flow sensor to activate the nonlinear numerical dissipation portion of a high order shock- 
capturing scheme ^j+i/v anc ^ K is a positive parameter that is less than or equal to one. The choice of 
the parameter k can be different for different flow types and is automatically chosen by using the local k 
described in [13]. However, in this study we set n = 1 to be used with the standard Durcos et al. flow sensor. 

It is noted that the nonlinear filter step described above should not be confused with the LES filtering 
operation. For previous studies on the performance of this filter scheme in DNS and LES simulations, 
see [14, 25, 26, 27, 28, 29, 13, 15]. This scheme has been validated for DNS of a 3D channel flow, 2D 
temporal and spatial evolving mixing layers, Richtmyer-Meshkov instability, 3D Taylor-Green vortex, 3D 
isotropic turbulence with shocklets, extreme condition flows, and 3D LES of temporal evolving mixing layers 
comparing with experimental data. 


2.3 Modifications of the LES Filtering Procedure for Flows with Shocks 

During LES computation using the filtered governing equations, there are two additional sources of inaccu- 
racy that may appear near the shock. The first one is connected with the numerical scheme used for solving 
the governing equations. Away from the shock the high-order central scheme is applied, introducing a neg- 
ligible amount of numerical dissipation. However, in the vicinity of the shock the shock-capturing scheme is 
activated, introducing numerical dissipation into the computed solutions. The amount of numerical dissipa- 
tion introduced by the shock capturing scheme depends on the particular problem parameters and may be 
higher than the turbulent dissipation modeled by the SGS terms. Hadjadj and Barmeo [7, 8] proposed to 
locally disable the SGS terms in order to obtain more accurate results. Denote this procedure as LES-Z: 


4 



Procedure LES-Z 

1. Employ flow sensor to detect shocks 

2. Obtain C s from Dynamic Germano Procedure except near the shock 

3. Set C s = 0 within grid stencil width of the scheme at shock 

Setting C s = 0 within the grid stencil width of the shock location (according to the computed flow sensor 
indicator) might appears to be too abrupt. An alternate procedure would be to set C s = 0 at the shock 
location and a smooth transition function C s such as tanh in the vicinity of the shock. Our studies show 
that setting C s = 0 for the whole WENO grid stencil width rather than a tanh smoothing function produces 
better results. 

The second additional source of inaccuracy of LES results in the vicinity of the shock comes from the 
fact that the explicit filtering operation in (10) - (13) is applied across the shock, causing inaccuracy of 
the results. In this case, as is pointed out in [2], a one-sided filtering operation should be used instead of a 
central form. Denote this procedure as LES-1S: 

Procedure LES-1S 

1. Employ flow sensor to detect shocks 

2. Employ Harten’s subcell resolution approach to detect exact shock location 

3. Obtain C s from Dynamic Germano Procedure except near the shock 

f. Apply one-sided filtering procedure to the left and right of the shock location 

It is noted that for the one-sided filtering formulation, one should use a discrete filter with the same filter 
width as in the standard central filter procedure. For example, with the trapezoidal filter (14) one can use 
the right-sided filter /+ = + \f i+ 1 , which according to formula (15) has the same A = h\ff> as filter 

(14). 

The third considered modification is based on Harten’s subcell resolution (SR) approach [6] combined 
with ENO reconstruction. Any shock-capturing method applied at the shock will result in smearing the 
shock. By estimating the exact shock location through SR and using ENO reconstruction one can sharpen 
the shock so that the numerical dissipation is minimized more locally. Denote this procedure as LES-SR: 

Procedure LES-SR 

1. Employ flow sensor to detect shocks 

2. Employ Harten’s SR approach to detect exact shock location 

3. Obtain C s from Dynamic Germano Procedure except near the shock 

j. Employ ENO reconstruction for the values of the points adjacent to the shock 

5. Apply one-sided filtering procedure to the left and right of the shock location 

3 3D Turbulence Across a Supersonic Shock Wave 

3.1 Problem Setup 

The 3D test case considered here concerns an initial turbulence disturbance at inflow boundary interacting 
with a stationary supersonic planar shock wave. The problem has been studied by previous investigators, 
mainly related to DNS computations, e.g. [30, 9, 8[. Here we choose the configuration considered in the 
DNS study of [9] . Figure 1 shows a schematic of the problem setup. The computational domain limits are 
—2<x< —2 + 47t, 0 < y < 2n and 0 < 2 < 2-7T. The grid is uniform in all directions with the spacing 
in x three times finer than in y and 2 (see [9] for explanation). We solve the filtered governing equations 
(1) - (3) in a non-dimensional form. The Yee et al. filter scheme with Ducros et al. flow sensor [24] is 
used for integration of the governing equations. The spatial base scheme is the 8th-order central differencing 
and the nonlinear filter scheme is the dissipative portion of the 7th-order WENO scheme. Since the initial 
data consists of a planar shock in the ^-direction, numerical dissipation should be mainly needed in the 
^-direction. In order to obtain more accurate results WENO dissipation is employed only in the x direction 
at the postprocessing stage of the Yee et al. filter. The inflow and outflow boundary conditions are applied 
in the streamwise direction and periodic boundary conditions are applied in the transverse directions. 

Inflow boundary condition. A fully developed turbulent inflow condition is applied using a turbulent 
database. This database is generated as follows. First, an initial isotropic turbulent field with the energy 
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spectrum E(k) ~ k 4 exp(— 2k 2 /kg) and microscale Reynolds number Re\ = pXu' rms / p = 140 is generated 
using the methodology described in [31]. Here the energy peak wavenumber ko = 4 is used. Next, the 
decay of this field in a periodic box is simulated for approximately three eddy turnover times r = A/ v! rms to 
ensure fully developed turbulence. After the decay the Reynolds number Re\ = 40 and the turbulent Mach 

number is M f = it'u' /cp = 0.16. Here Co is the mean speed of sound. The generated isotropic turbulence 
is introduced at the inflow boundary with constant mean velocity up. We consider two cases with mean flow 
Mach number M = 1.5 and M = 3.0. In order to compare the DNS results we use the inflow database from 

[91- 


Outflow boundary condition. In order to avoid acoustic reflections of subsonic flow from the outflow 
boundary a non-reflective sponge layer is employed on the region near the outflow. The length of this layer is 
x max ~ Xsp = 7r . The sponge layer is implemented by introducing a following source term into the equations 


(1) - (3) 


ft = — 


kpiip 
27 r 


X %sp 
% max 


(/ < / >yz)t 


(19) 


where f = p,piLi,pE and {-) yz denotes averaging in the y- and z- directions. 

The outflow pressure p ^ is chosen such that the mean shock location is stationary. For laminar flow 
Rankine-Hugoniot conditions give 


Poo = 1 2 7 ~ (u 0 - U s ) 2 

Po 7 + 1 L Cp 


(20) 


where po is the inflow mean pressure, Up is the mean inflow velocity, Cp is the mean inflow speed of sound 
and U s is the shock velocity. As the inflow condition is turbulent, the Rankine-Hugoniot conditions are valid 
only instantaneously but not in average. After an initial guess based on (20) the outflow pressure is refined 
by an iteration procedure, integrating the governing equations on a coarse grid and updating the pressure 
according to the formula 

P'oo = Poo + 4U s p u Up/ (7+1)- (21) 

See [9] for more details. 

Gathering statistics. The simulation statistics for a given function / are obtained by averaging in time 
and in homogeneous directions: 


< f > H (x) 


1 

LyL z /S.t 




f(r,t)dzdydt, 


(22) 


where L y and L z are domain sizes in transverse directions. The averaging is performed over a time span 
At « 100/(/coUo)- Convergence is confirmed by comparing the results with statistics obtained over time span 
At/2. The statistics are gathered after the transition period has passed. Transient time tp is estimated as 
to L x /up, where L x is the domain size in the streamwise direction. The correct choice of transient time 
is confirmed by comparing with the statistics obtained starting from time to/2. 


3.2 DNS Comparison 

After a grid refinement study the DNS results as a reference solution for LES comparison are obtained on 
a grid with 1553 x 256 2 points. The instantaneous streamwise and transverse velocity fields are shown in 
Fig. 2. During the computation over a long time evolution, the shock slightly moved upstream. For this 
turbulent Mach number the shock is wrinkled due to turbulent inflow. As shown in previous studies [30, 9], 
the shock may break at higher turbulent Mach numbers M t . The turbulence is compressed by the shock, 
and immediately behind the shock it is anisotropic. The comparison of turbulent statistics for streamwise 
and transverse components of the vorticity and Reynolds stress shows that the turbulence becomes isotropic 
again downstream of x sa 3. Downstream of x = 8 turbulence is essentially damped with the sponge source 
term. 

Figures 3-6 represent the DNS statistics for different cases. The plots are presented in the reference 
frame of the average shock location. 

The comparison of the DNS statistics obtained in this work using the ADPDIS3D [5] code with the 
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Post-shock Zone Sponge Layer 



Figure 1: Turbulence across a shock wave problem setup. 



Figure 2: Instantaneous velocity field u x (top) and u y (bottom) obtained with DNS on grid of 1553 x 256 2 
points. Slice z = const. 
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data obtained from digitizing the solution [9] are shown in Fig. 3. The results are in good agreement. 
The grid resolution in the vicinity of the shock is the same. In [9] grid clustering near the shock has been 
employed, resulting in a smaller grid size, 694 x 256 2 . The work [9] employs the HYBRID code, which also 
uses a Ducros et al. flow sensor for shock detection. However, in [9] WENO dissipation has been introduced 
at every Runge-Kutta stage, whereas the Yee et al scheme allows decreasing the computational cost by 
employing WENO dissipation only after the full Runge-Kutta step. 

The results of the grid refinement study for cases M = 1.5 and M = 3 on levels of grid refinement are 
shown in Figs. 4 and 5. For the case M = 1.5 the results on 777 x 128 2 are fairly close to the DNS results 
obtained with 1553 x 256 2 grid points, especially for the Reynolds stress statistics. For the case M = 3 the 
results also converge, but the convergence is slower. 

The comparison between the filtered DNS data and the results obtained with different numerical schemes 
on the coarse grid with 389 x 64 2 points is shown in Fig. 6. Here we compare the results obtained by WENO 
of 5 th order (WEN05) and WENO of 7 th order (WEN07) with the results obtained by the filter counterpart 
of WEN07 with Ducros et al. split (WEN07fi+split). For DNS on a coarse grid a too low-dissipative scheme 
will produce inaccurate results. To obtain the results by the filter scheme shown in Fig. 6, we use WENO 
dissipation in all directions. It is noted that for LES simulation it is more important to use a low-dissipative 
scheme, since the dissipation is already introduced by the turbulence model. 

3.3 LES Comparison 

In this section we compare the results obtained by the Germano model using different filtering procedures 
(LES with standard filtering procedure, LES-Z, and LES-1S) with the DNS solution filtered to the size of 
the LES grid. Results obtained by LES-SR are not shown, since that method is still under development. 

Comparison of the methods on a grid with 389 x 64 2 points for the case M = 1.5 is shown in Figs. 8 and 
9. For this case the results obtained using LES-Z and LES-1S are closer to the filtered DNS than standard 
LES. For certain variables, LES-1S performs slightly better than LES-Z. However, the difference between 
LES-1S and LES-Z may be more significant for other flow conditions, e.g., higher Mach and turbulent Mach 
numbers. In the case when the shock-capturing scheme dissipation might be larger than the SGS dissipation, 
it is conjectured that LES-Z would be more accurate than the other considered approach. But for the cases 
where the SGS dissipation is larger than the numerical scheme dissipation, LES-1S might obtain more 
accurate results than LES-Z. 

The results for the same case on a two-times finer grid ( 777 x 128 2 ) are shown in Fig. 10. All results 
obtained on this grid are quite close to each other, though LES without any modification of the filtering 
procedure is slightly closer to the filtered DNS data. One can assume that the shock is resolved enough at 
this level of grid refinement so that there is no need for additional shock treatment. 

A similar comparison of the methods on the grid with 389 x 64 2 points for the higher Mach number case 
( M = 3) is shown in Fig. 11. The behavior of the methods is similar to the M = 1.5 case, but LES-Z here 
performs slightly better than LES-1S, which can be seen especially in the plot of the y vorticity component. 

It is remarked that as we employ implicitly filtered LES equations, the width of implicit filter is not well 
defined. One should take this fact into account when using filtered DNS data as a reference solution for 
verification of the LES model. Figure 12 shows the differences in filtered DNS data when using two different 
filters: 

• the 5-point top-hat filter with coefficients (1/8, 1/4, 1/4, 1/4, 1/8) and width A « 4.2 h according to the 
estimate (15), 

• the 7-point trapezoidal filter with coefficients (1/22, 1/11, 2/11, 4/11, 2/11, 1/11, 1/22) and width A « 
4.8 h. 

The top-hat filter has been used to obtain DNS filtered data plotted on previous figures. For more precise 
comparison of LES filtering modifications, explicit LES will be considered in the future. 


4 Conclusions 

The DNS results obtained by high order nonlinear filter schemes compare well with the reference solution 
[9]. In general, the employment of the Yee & Sjogreen filter schemes requires less computational cost than 



standard Hybrid schemes. Our LES study confirms the results found by previous authors that the dynamic 
Germano LES model with a standard filtering procedure may loose accuracy due to a strong shock. Two 
modifications of the LES filtering procedure (LES-Z and LES-1S) which are designed for improving the 
accuracy of the standard method have been considered. For this particular shock-turbulence interaction 
test case both modifications of the filtering procedure show similar results which are more accurate than the 
results obtained using the standard LES filtering procedure. However, the turbulent Mach number M t = 0.16 
for the test case considered here is quite low, and the SGS dissipation may be not high enough in comparison 
with numerical dissipation of the shock-capturing scheme. Different behavior of considered procedures is 
expected for high turbulent Mach numbers, which is forthcoming. In addition, a systematic assessment 
employing LES-SR and LES-1S will be reported in a forthcoming report. Also, for better estimation of the 
simulations accuracy the performance of explicitly filtered LES will be considered. 
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— ADPDIS3D 
+ HYBRID 


Figure 3: Comparison of DNS statistics obtained in this work employing the ADPDIS3D code with data 
obtained from digitizing the solution [9] employing the HYBRID code. Top row: Reynolds Stress components 
x (left) and y (right). Bottom row: vorticity components x (left) and y (right). 



Figure 4: Grid refinement study for the case M = 1.5 on three grid levels: 389 x 64 2 , 777 x 128 2 and 
1553 x 256 2 . Top row: Reynolds Stress components x (left) and y (right). Bottom row: vorticity components 
x (left) and y (right). No LES modeling. 
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Figure 5: Grid refinement study for the case M = 3 on three grid levels: 389x 64 2 , 777 x 128 2 and 1553 x 256 2 . 
Top row: Reynolds Stress components x (left) and y (right). Bottom row: vorticity components x (left) and 
y (right). No LES modeling. 





Filtered DNS 

WEN07fi+split 

— WEN07 
WEN05 


Figure 6: Comparison to filtered DNS data of the statistics obtained by the different numerical schemes 
(WEN05, WEN07 and WEN07fi+split) on grid 389 x 64 2 , M = 1.5. Top row: Reynolds Stress components 
x (left) and y (right). Bottom row: vorticity components x (left) and y (right). No LES modeling. 
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Figure 7: Comparison to filtered DNS data of the statistics obtained by LES using different numerical 
schemes (WEN05, WEN07 and WEN07fi+split) on grid 389 x 64 2 , M = 1.5. Top row: Reynolds Stress 
components x (left) and y (right). Bottom row: vorticity components x (left) and y (right). 





Filtered DNS 

No LES model 
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— LES-1S 


Figure 8: Comparison to filtered DNS data of the statistics obtained by LES with different filtering procedures 
(standard LES, LES-Z and LES-1S) on grid 389 x 64 2 , M = 1.5. Top row: Reynolds Stress components x 
(left) and y (right). Bottom row: vorticity components x (left) and y (right). 
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Figure 9: Same as Figure 8, zoom in the vicinity of the shock 





Filtered DNS 

No LES model 
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Figure 10: Comparison to filtered DNS data of the statistics obtained by LES with different filtering pro- 
cedures (standard LES, LES-Z and LES-1S) on grid 777 x 128 2 , M = 1.5. Top row: Reynolds Stress 
components x (left) and y (right). Bottom row: vorticity components x (left) and y (right). 
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Figure 11: Comparison to filtered DNS data of the statistics obtained by LES with different filtering proce- 
dures (standard LES, LES-Z and LES-1S) on grid 389 x 64 2 , M = 3. Top row: Reynolds Stress components 
x (left) and y (right). Bottom row: vorticity components x (left) and y (right). 





DNS: Top-hat 5 pt filter 

— DNS: Trapezoidal 7pt filter 
LES 


Figure 12: Comparison to filtered DNS data using different filter operator (Trapezoidal and Top-hat). Grid 
389 x 64 2 , M = 3. Top row: Reynolds Stress components x (left) and y (right). Bottom row: vorticity 
components x (left) and y (right). 
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